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Here we report on the design, construction and characterization of a time domain Balanced Homodyne Detection 
apparatus operating in the regime of high electronic noise. The experimental setup has been commissioned by 
studying coherent ultrashort (~ 80 fs) light pulses using Pattern-function Quantum Tomography. Furthermore, 
we proved that the electronic noise can be treated as a detector inefficiency showing the effectiveness of this 
approach by measuring light coherent states with different energies. 



In Quantum Mechanics the measurement of the state 
of a physical system implies to estimate the expectation 
value of any observable of the system. For a generic 
operator the expectation value is obtained starting from 
the measurement of an ensemble of observables defined 
"quorum''^. In particular, for the electromagnetic field 
the quorum is a continuum set of the field quadra- 
ture^. 

Balanced Homodyne Detection (BHD) is an experimen- 
tal method for measuring a discrete set of the field 
quadratures for a single mode cw fieldStH. Such mea- 
surements can be analyzed by the Pattern-function 
Quantum Tomography (PFQT) technique in order to 
characterize the quantum state of the optical mode^. 
Instead, for pulsed radiation fields a BHD operating 
in the time domain must be usedl^^^. In this case the 
optical pulses are equally prepared quantum states and 
the set of measurements fed to PFQT comprises single 
pulse quadrature measurements. Although this is a 
standard method in quantum information, its effective- 
ness under high noise condition has to be proven. 
Here we report the characterization of a time domain 
BHD apparatus operating in a regime of high electronic 
noise (EN). We treat the low shot-to-electronic-noise 
ratio (w 0.8 dB) as a detector inefficiencytlS]. In order 
to prove the applicability of this treatment to high EN 
homodyne detection, the PFQT technique has been ex- 
ploited. The expectation value of the number operator 
of pulsed coherent states with different energies has 
been found consistent with the measured mean number 
of photons per pulse. Furthermore, the analyzed states 
are confirmed to be Gaussian and minimum uncertainty 
states. 



Results 

In BHD the signal, i. e. the field under investigation, is 
mixed with a strong coherent reference field, the local 
oscillator (LO), by a 50/50 beam splitter. The outputs 
are collected by two photodiodes and the difference 
photocurrent (homodyne photocurrent) is measured. It 
can be proven that, when the is significantly more 
intense than the signal, the homodyne photocurrent is 
proportional to the signal field quadrature!^. Given a 
and a} the field operators associated with the optical 
mode of the signal, the quadrature operator is defined 
as: 



where $ is the relative phase between the signal and the 
LO. The continuum set of quadratures with $ e [0, 2tt] 
constitutes a quorum characterizing the signal mode. 
Noteworthy the concept of light mode, usually associ- 
ated with a monochromatic field, needs to be general- 
ized for pulsed light. A thorough theoretical treatment 
for BHD in pulsed regime is presented herewith in the 
Methods section. 

The opto-mechanical scheme of our experimental ap- 
paratus is shown in Fig. ^ The laser source is a 
mode-locked Ti:Sapphire oscillator with 80 MHz repeti- 
tion rate. A beam splitter divides the incoming beam 
in two parts which then interfere in a second beam split- 
ter (NPBS in Fig. [TJ. The outputs are detected and 
subtracted by a differential photodetector ( Thorlabs 
PDB430A). The latter is made up of two Si /FIN phcn 
todiodes with nominal quantum efficiency of 0.77. We 
verified that the response of the photodiodes is linear 
up to 0.6 mW LO power, while a non-linear response 
sets in at higher powers. The detector subtraction ef- 
ficiency is quantified by the common mode rejection 
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Figure 1: Scheme of the opto- mechanical setup. 



ratio (CMRR) , defined as the ratio between the detector 
output power when both photodiodes are illuminated 
and the power when one of the two is screened. The 
measured CMRR for our apparatus is larger than 36 
dB. 

The homodyne photocurrent is recorded by a digital 
oscilloscope ( Tektronix TDS3000B) with a bandwidth 
of 500 MHz and a sampling rate of 5 GSamples/s. The 
digitized output is numerically integrated over time in- 
tervals corresponding to the duration of the pulse. Each 
integral is associated with a single quadrature measure- 
ment. We checked the independence of each quadrature 
measurement by performing a correlation test between 
two subsequent pulses. The Correlation Coefficients 
(CCs)'Sl between pulse n and n + m (m = 0, 1...9) are 
shown in Fig. |2](a). In the inset a parametric plot of 
2 X IC* subsequent pulses is shown, where the pulse 
is plotted against the pulse n. The weak correlation 
in this plot demonstrates that there is no significant 
impact on the measured integral of the pulse n + 1 
from that of the pulse n. Each CC in Fig. [2] (a) is 
calculated using the data obtained from 200 difference 
pulses. The CC between pulse n and pulse n -I- 1 is 
0.30 ± 0.08; this value is slightly higher with respect to 
other reported BHDs, i. e. smaller than O.fl^^HIIl. We 
can conclude that there is a small correlation between 
adjacent pulses, since the CC has approximately the 
same value for m = 2, ...9. This could be attributed to 
noise at a frequency lower than the repetition rate of 
the pulses. 

In shot noise regime the homodyne detector (HD) noise 
variance is expected to change linearly in the LO power 
with a constant offset representing the electronic noisel^^. 
We calculated the HD noise variance of 8 x 10'^ pulses 
for different values of the LO power. The result is 
shown in Fig. |2] (b): the HD noise variance contains 
an electronic background and a linear contribution up 
to 0.6 mW LO power; for higher powers the non-linear 
effects due to the photodiodes start. We work at 0.6 
mW LO oscillator power in order to have the maximum 
shot-to-electronic noise ratio (S* — 1.2) achievable in 
the linear regime. The measurements reported in Fig. 
[2] are obtained with the signal beam blocked, i.e. with 
the signal in the vacuum state. 

When the signal beam is not blocked, it is attenuated 
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Figure 2: Detector characterization measurements, (a) 
Correlation coefficient between pulse n and n + m {LO 
power: 0.6 mW). Inset: parametric plot of 2 x 10^ 
subsequent pulses where the pulse n+1 is plotted against 
the pulse n. (b) HD noise variance versus LO power; 
shot noise contribution (dashed curve); electronic noise 
background (dashed-dotted curve). 



with respect to the LO by means of a neutral density 
filter {Fi in Fig. [ij. The phase difference $ between the 
two arms of the interferometer can be modulated in the 
[0, 2tt] range using a piezoelectric translator in the LO 
arm (PZT in Fig. [l]). The homodyne photocurrent is 
acquired by the digital oscilloscope and each difference 
pulse is integrated. The obtained value Ni, correspond- 
ing to a certain piezo position pi , is proportional to a 
quadrature measurement Xt^: 



N, 



(2) 



The proportionality factor 7 is obtained using the vac- 
uum state as a reference. This is allowed since the 
quadrature variance for the vacuum state is ^^[x^Jio) = 
1/2 for any phase value. Defining 7 as: 



(3) 



the vacuum homodyne measurement has been used to 
extract the Nio values. (N^q) is the variance of 8 x 10^ 
experimental data of the vacuum. The homodyne data 
have been calibrated dividing Ni by 7. Note that 7 will 
give a misleading estimation of the actual mean photon 
number. Therefore, a correction taking into account 
the EN effect is necessary, as discussed in the following. 
The calibrated measurements for three optical coherent 
states with different energies are reported in Fig. |3] 
The optical density CD of the filter Fi used for each 
measurement is given in the caption. The homodyne 
traces are composed of M = 8 x lO'' experimental data: 
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Figure 3: Homodyne traces of three optical coherent 
states. Each trace is acquired using a different optical 
density OD of the filter Fi. (a) OD = 4.5; (b) OD = 7; 
(c) OD = 8. In the inset of (a) we show the interference 
figure obtained by measuring the mean value of 4 inte- 
grated pulses versus the piezo position. The homodyne 
traces are measured in the piezo range corresponding 
to the central optical cycle of the interference figure. 

{(xi; $i), (x2; $2), -, {xk;<^k), {xm] <^m)}- 
However, the effect of the EN, which adds a random 
quantity to each field quadrature measurement, must 
be accounted for. Appel at al. demonstratel^^ that this 
noise can be treated as an optical loss channel with 
equivalent transmission efficiency, 



eq 



(4) 



where S is the ratio between shot and electronic noise at 
the chosen LO power. In particular, the vacuum state 
measurement is distorted by the EN and consequently 
the rescaling factor 7 must be changed. Following the 
Appel's treatment's] and considering equation the 
following expression is obtained. 
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where 7' is the conversion factor in case of absence of 
electronic noise. Thus, in order to take into account 
the EN, the homodyne data have been normalized to 7'. 



These data have been used to characterize the quantum 
state of the signal field via PFQT^^. 



Discussions 

The PFQT technique can be summarized as follows. 
For a generic observable O of the optical mode it is 
possible to find a pattern function Rq{x; <&) for which 
the statistical average over the experimental data. 



(6) 



k = l 

gives an estimation of the expectation value (O) with 
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This technique is used to estimate the expectation value 
of the number operator {a^a), which can be evaluated 
in an independent way. 

The pattern function of the number operator is: 
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where the factor 77 is the quantum efficiency of the 
homodyne detection apparatus. In order to take into 
account the EN loss effect, we put 77 = rj^q given in 
equation Q. Finally we estimate the mean photon 
number as: 
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(9) 



The error on this estimation is calculated using equa- 
tion Q. The results for the three different homodyne 
traces are reported in Table [T] The expectation values 
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4.5 


596 ±2 
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1.6 ±0.1 
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0.6 ±0.1 



Table 1: Expectation values of the number operator of 
three optical coherent states for three different optical 
densities ODs. The values are obtained via PFQT from 
the homodyne measurements reported in Fig. |3] taking 
into account the equivalent quantum efficiency r]eq as 
in equation Q. 

calculated using the PFQT are consistent, in the order 
of magnitude, with the mean number of photons per 
pulse obtained from the nominal OD of the filter. In 
this way we verify that the treatment of low EN'^ is 
valid also for high EN detection systems. 
Now we focus on the optical state with the lowest en- 
ergy [Fig. [3] (c)]. The analysis for the other states is 
analogous. 
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First we verify the Gaussianity of the experimental data 
through the evaluation of the kurtosis excess k, which 
is zero for Gaussian distributionsli^l. Considering 10'^ 
experimental data in the first phase range of the ho- 
modyne trace in Fig. |3] (c) and we obtain a kurtosis 
excess k = —0.008. Since all the phase ranges lead to 
a K ~ 0, this proves the Gaussianity of the homodyne 
probability distribution function. 

Then we use the PFQT to evaluate the expectation 
values of the the two quadratures associated with 
the position and momentum operators (aij^^o) = 9; 
a;(*=7r/2) = p) s-iid of their variances (cr^[(7]; (T^[p]). The 
results, shown in Table [2j demonstrate that the ana- 
lyzed state is a minimum uncertainty state. 
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0.36 ±0.01 


0.53 ±0.05 


(P) 




0.00 ±0.01 


0.53 ±0.05 



Table 2: Expectation values of the position and mo- 
mentum operators and their variances. The values are 
obtained via PFQT from the homodyne measurements 
reported in Fig. [3](c). 

Finally we evaluate the covariance matrix C, that fully 
characterizes Gaussian states and directly depends on 
the field quadratures'^. The elements of the covariance 
matrix are defined as: 
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Cfcj = 7^{{Rk,Rj}) - (Rj) (Rk) 



(10) 



where R — {q,p) is the first-moments vector of the opti- 
cal state and {A, B} — AB+BA is the anti-commutator 
between two operators From the PFQT we obtain 
the following estimations: 



C 



0.53 ±0.5 
-0.02 ±0.03 



-0.02 ±0.03 
0.53 ±0.05 



(11) 



This is consistent with what is expected for a coherent 
state. The other two homodyne measurements, 
reported in Fig. [3] (a) and (b), lead to similar results. 

In conclusion, we reported the characterization of 
a time domain BHD apparatus operating in the high 
electronic noise regime. The apparatus has been char- 
acterized via the PFQT of optical coherent states with 
different mean number of photons per pulse. We treated 
the electronic noise trough the approach used by Appel 
at al.ESl and we showed that such a treatment is valid 
for high EN homodyne detectors. This result is of rele- 
vance in the perspective of using BHD and PFQT in 
experiments performed in high noise environments such 
as time domain pump-probe spectroscopy. In particular, 
ultrashort light pulses could be used to prepare coherent 
states into matter and BHD to have a full character- 
ization of the light pulses which have interacted with 
transient matter states. 



Methods 

BHD theoretical treatment in the pulsed regime. 

BHD in a pulsed regime requires a formal generalization 
of its theoretical description with respect to the one- 
mode regime; indeed, the LO and the signal at the 
beam splitter are not monochromatic. Classically, the 
electric field of a pulsed laser beam is in a mode-locked 
superposition of amplitudes: 
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'S?i{t)^LUit + ipi , (12) 



where the phases ^i{t) are mode- locked by the condi- 
tion ifi — lipo- Each frequency contributes to the 
field with the amplitude a; — \ai \ e'*' and the number 
of contributing frequencies depends on the shape of 
the pulse. Quantized pulsed laser light is described 
by associating to each monochromatic component a 
coherent state |a/), that is an eigenstate of the annihi- 
lation operator a/ of photons in the mode of frequency 
uji, ai \ai) = ai \ai) , and to the entire pulse the tensor 
product 



M 



(13) 



l=-M 



where a is the vector whose components are the am- 
plitudes ai. By means of the creation and annihilation 
operators a and dj each monochromatic coherent state 
reads 



\ai)^D{ai)\0) , i?(aO-e"''' , 



(14) 



where |0) is the vacuum state and D{a) is the so-called 
displacement operator. 

Since the creation and annihilation operators pertaining 
to different modes commute, the pulsed coherent state 
(13) can conveniently be recast as: 



\a) = D(a) |0) , D{a) = e 



A'{a)-A{a) 



(15) 



in terms of a multi-mode displacement operator D{a) 
where 



A\a) — ai a| , A{a) — a;*d; 



(16) 



The reason for labeling the pulsed coherent state by \a) 
can now be easily understood. The state of one photon 
of frequency ui is given by d\ |0) = jl;), while a generic 
non-monochromatic superposition of frequencies uii with 
amplitudes a; corresponds to the state jig) — ai \ 
that results from 



it(a)|0)=5]a,aj|0) 



E 



II. 



(17) 



Therefore, A^{a) is the creator operator of a single- 
photon in the (not normalized) superposition state |la), 
while A{a) destroys a photon in the same state; thus, 
the quantum state of the pulsed laser is a coherent 
state associated not with a single amplitude ai, but 
with the vector a of all the amplitudes contributing 
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to the pulse: in other words, we have a Poissonian 
distribution not with respect to the number of photons 
in a monochromatic wave, but to the number of photons 
in the superposition |la). 
The normalized operators 



A = 



it 



A{a) 



(18) 



where |ap = (la|la), satisfy the canonical commuta- 
tion relations 



i,it 
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^a*aj [ai,a]] = 1 . (19) 



For this reason, in the main text of the letter we have 
limited the discussion to the simpler case of a single 
mode BHD, with the proviso that whenever a quadra- 
ture is used, it actually refers to its expression in the 
pulsed regime. In particular, the measured data can be 
used to reconstruct the expectation values of all signal 
observables that can be expressed as functions of the 
operators A and . For instance, the mean value of 
the photo-current operator second moment with respect 
to the LO pulsed coherent state \z) {z\ is 
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(24) 



By means of them we can now extend the state- 
tomography techniques to the case in which the signal 
and the LO are pulsed. In the monochromatic case, 
when the signal mode a interferes with the LO mode 
b on the beam splitter, the photo-current operator is 
J = at6-h&tai7^ 

If more frequencies uji are present both in the LO and in 
the signal, each one of the corresponding mode operators 
will be subjected to the beam splitting transformation 
and the detectors will ideally register photons of all 
involved frequencies. Then, the photo-current operator 

becomes I = J2i ^l^i ~^ 

If the LO incoming field is in a pulsed coherent state 
1^) = e'^*'^-'"^*^^^ |0), with generalized creation and 
annihilation operators B^{z) = J2i ^^"^ ^(^) — 
J2i^i^h the phase difference $ between the LO and 
the signal is changed by the action of the piezoelectric 
translator placed in the LO arm on all the LO modes: 

^ b] e-** . 



bi e'* 



(20) 



The photo-current operator which is measured by the 
pulsed honiodyne setup is thus given by 



U^^(a\bi e'* + aib\e 



(21) 



Let ps denote the quantum state (density matrix) of 
the signal field and by \z) {z\ the projector onto the 
(coherent) state of the incoming pulse. Then, using 
that bi\z) — zi \z) and the expressions in (16), the ex- 
pectation value /$ — Tr (^ps i8) \z) 
current is calculated as follows: 
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Using ( 18 \ and comparing these results in the standard 



treatment of BHD for the single-mode case, one realizes 
that the pulsed BHD measures a quantity propor- 
tional to a quadrature which generalizes that in ([T]): 



/$ = \/2|z|Tr 



Xif, 
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(23) 



and differs from the pulsed quadrature second moment 
X|, by the number of photons in the pulsed signal which 
is to be taken much smaller than the intensity |zp of 
the pulsed LO. A similar suppression by occurs for 
the correction terms appearing in higher moments so 
that the distribution of the outcomes of the homodyne 
photocurrent is equal to that of the corresponding field 
quadratures. 

Notice that, in the case of only one frequency mode, 
the above setting reduces to the standard homodyne 
tomographic techniques. Therefore, in conclusion, we 
have generalized the treatment of a monochromatic 
BHD to the case of pulsed input light. It explicitly 
takes into account the fact that the state of a pulsed 
laser is the tensor product of the coherent states of each 
contributing frequency. 
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